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Abstract 

The Jacobian matrix of a dynamic system and its principal mi- 
nors play a prominent role in the study of qualitative dynamics and 
bifurcation analysis. When interpreting the Jacobian as an adjacency 
matrix of an interaction graph, its principal minors correspond to sets 
of disjoint cycles in this graph and conditions for various dynamic 
behaviors can be inferred from its cycle structure. For chemical re- 
action systems, more fine-grained analyses are possible by studying a 
bipartite species-reaction graph. Several results on injectivity, multi- 
stationarity, and bifurcations of a chemical reaction system have been 
derived by using various definitions of such bipartite graph. Here, we 
present a new definition of the species-reaction graph that more di- 
rectly connects the cycle structure with determinant expansion terms, 
principal minors, and the coefficients of the characteristic polynomial 
and encompasses previous graph constructions as special cases. This 
graph has a direct relation to the interaction graph, and properties of 
cycles and sub-graphs can be translated in both directions. A sim- 
ple equivalence relation enables to decompose determinant expansions 
more directly and allows simpler and more direct proofs of previous 
results. 

1 Introduction 

A major problem in the analysis of chemical reaction systems is the fact 
that parameters such as kinetic rate constants for are inherently difficult 
to obtain from experimental data, and in-vitro data might not be valid 
for in-vivo systems. Moreover, the exact algebraic form of chemical rate 
laws is often difficult to determine. In contrast to more general dynamic 
systems, however, the dynamics of a chemical reaction system is addition- 
ally restricted by the stoichiometry and topology of the underlying reaction 
network, and chemical rate laws are typically monotone in the concentra- 
tions. The underlying reaction network also induces algebraic dependencies 



of state variables and parameters. Several approaches have been proposed 
that exploit these additional constraints and allow to determine — without 
knowledge of parameter values and minimal conditions on the rate laws — if 
a system is capable of certain qualitative dynamics such as oscillations and 
multiple equilibria, and to establish stability of equilibria. These methods 
fall broadly into three categories: 

(i) methods based on the particular algebraic structure of reaction sys- 
tems include the deficiency-based Chemical Reaction Network Theory [15} 
[TOl [9], the Stoichiometric Network Analysis [4], and Biochemical Systems 
Theory [El EH [20] . 

(ii) methods based on the cycles of a signed interaction graph derived 
from the Jacobian: e.g., a positive cycle is necessary for multistationar- 
ity, whereas a negative cycle is necessary for oscillations [HJ [2H [231 US]- 
Other criteria have also been investigated [S] . The absence of negative undi- 
rected cycles implies monotonicity with respect to an orthant cone, ruling 
out chaotic and oscillatory dynamcis [22l [211 E! ! 

(iii) methods based on a bipartite graph with vertices for both species and 
reactions been proposed recently to study criteria for the existence multiple 
equilibria or oscillations [31 H7J El [7] and for establishing monotonicity [22], 
[1]. Exploiting the relation of principal minors of a matrix and cycles in 
its associated graph to coefficients in the characteristic polynomial, several 
criteria were also formulated to infer the possibility of saddle-point and 
Hopf bifurcations in a system |17| . Definitions of such graphs in the context 
of qualitative dynamics include undirected [7J [3] and directed [IT] species- 
reaction (SR) graphs as well as graphs with special edge-types For 
various proofs, oriented versions of an undirected SR-graph also need to be 
considered in [7J. 

In this work, our goal is to elucidate how the various interaction- and 
bipartite-graph-based methods are related and to propose a new version 
of the species-reaction graph, named directed species-reaction graph (DSR- 
graph), that immediately shows its relation to the Jacobian matrix and 
gives more direct justification and proofs of several theorems. Using this 
graph, we present refined criteria for sign-definiteness of the determinant of 
the Jacobian and its principal minors and answer a question raised in [T7J 
of the equivalence of criteria developed independently in [17] and [7J. We 
propose a simple equivalence relation on the DSR-graph that allows to study 
determinants by studying each of its equivalence classes. We also develop 
and emphasize the relation between determinant expansions, the Jacobian 
matrix, and the DSR-graph, making use of long-known relations [12] . 
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2 Chemical Reaction Systems 



We consider a biochemical reaction network with n species S\, . . . , S n and 
r reactions R\ , . . . , R r . A reaction Rj is formally given by 

n n 
i=l i=l 

where yi,j,y\j G K>o are the molecularities of the substrates and products, 
respectively. The matrix N G M nxr with entries Nij = y[ j—yi,j is called the 
stoichiometric matrix, its jth column vector (y[ ■ — yij, . . . , y' n • — y n ,j) T £ 
M n is the stoichiometry of reaction We require that all reactions are 
irreversible; this can always be achieved by splitting reversible reactions 
into a forward- and a backward reaction. We also require that each species 
occurs on at most one side in each reaction, i.e., we exclude reactions such 
as A + B — > 2A. The signs of all Ni j are then uniquely defined such that 
Nij > (resp. iVy < 0) if species Si is produced (resp. consumed) by 
reaction Rj, and Nij = if it does not participate in the reaction. 

We denote by Xi(t) G M>o the concentration of species S{ at time t, and 
collect these concentrations in the state vector x = x(i) = (xi(t), . . . , x n (t)) T , 
adopting the usual convention to drop the explicit dependence on t. Each 
reaction Rj is assigned a rate law 

Vj : M> — > M>o , 

describing the velocity of the reaction as a function of the current state x of 
the system. A typical rate law is the mass-action law Uj(x) = kj-x\ 1,j ■ ■ ■ x^' j 
with some reaction-specific rate constant kj > 0. The rate laws are collected 
into the rate vector v(x) = (vi(x), . . . , tv(x)) T , leading to the system of non- 
linear differential equations 

|x = f (x) = N • v(x) , (1) 

which govern the evolution of species-concentrations over time. 

A special class of rate laws is given by the non-autocatalytic (NAC) laws, 
for which dvj/dxi > when specie Si is a substrate of reaction Rj and all 
concentrations are positive, and dvj/dxi = if Si is not a substrate of 
Rj. These conditions are quite natural: they require that an increase in a 
reactant's concentration cannot decrease the reaction rate and that species 
that are not reactants of a reaction do not directly influence its rate. The 
class NAC contains both mass-action and Michaelis-Menten rate laws, as 
well as many others. Very similar conditions for reversible reactions were 
given in [3]. 

To model mass transport over the boundaries of the system, we also 
allow reactions of the form Rj : — > S{ with constant rate Vj(x) = kj > 
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to represent production or inflow of Si, and of the form Rj : Si — > with 
rate fj(x) = kj ■ X{ to represent degradation or outflow of 5$. 

A concentration x* G M™ is called an equilibrium or a steady state if 

f(x*) = ; 

it is called positive, if x* > for all i = l,...,n, which we abbreviate 
x > 0. With initial condition xo = x(0), all solutions of system ([!]) are 
confined to the afhne subspace xo + ImgN of dimension rank(N). A system 
is said to have multiple positive equilibria or to be multistationary, if there 
are x* / x** > such that 

f(x*) = 
f (x**) = 
x** - x* G Img N . 

Without loss of generality, we may assume that the system has been reduced 
such that N has full rank; then, ImgN coincides with the state space and 
the third condition for multistationarity is trivially fulfilled. 

Example 2.1. We consider the following reaction network 

R x : aA-> bB 
R 2 : cB -> dA 
R 3 : 
Ra : ->■ A 
R 5 : B^ 
R 6 : -> B . 

This network has n = 2 species A, B with concentrations x = (xa,xb) T 
and r = 6 reactions R%, . . . , Rq, where reactions R$ and R$ describe the 
degradation or outflow of species A, B, respectively, and reactions R4 and 
Rq their production or inflow into the system. Reactions R± and R2 are the 
internal or true reactions. The stoichiometry of reaction R\ is determined by 
2/1,1 = a : 2/2,1 = 0, 2/1 1 = 0, 1 = b and is thus (—a, b) T . The stoichiometric 
matrix is 

_ (-a d -1 1 0\ 
V& -c 0-1 l) 

with rank(N) = 2. The dynamics is given by the two differential equations 

(d/dt)xA = — a vi(x) + dt>2(x) — «3(x) + t>4(x) 
(d/dt)xB = -cu2(x) + 6^i(x) - tj 5 (x) + t> 6 (x) . 

The production rates U4(x) = £4 and t>6(x) = /c6 are constant. 
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3 Qualitative Dynamics and the Interaction Graph 
and 



Several criteria for qualitative dynamics such as oscillations or multistation- 
arity can be established via properties of the Jacobian matrix of a reaction 
system. The Jacobian can be interpreted in terms of a graph and its prin- 
cipal minors then correspond to certain cycle structures in this graph. Of 
particular interest are conditions such that these cycles are sign-definite, 
that is, do not depend on the concentrations x at which the Jacobian is 
evaluated. 

The Jacobian matrix of a reaction system ([!]) is the function 

,.(«) _ N .(*i 

\° x j / \<i,j<n \° x j / l<i<r,l<j<n 

Evaluated at Xq, the resulting matrix J(xo) allows to approximate the sys- 
tem's dynamics in the vicinity of Xq. 

A tupel a = (ai, . . . , a{) C {1, . . . , n} with a\ < ■ ■ ■ < a\ is called a 
multi-index of {l,...,n} of size \a\ := I. Let A £ R nxm and let a and 
f3 be two multi-indices, not necessarily of the same size, of {1, . . . ,n} and 
{1, . . . , m}, respectively. The matrix 

Aq.^ : = (a ) ; j j)ig Q , i jg i g 

is the | a | x \/3\ sub-matrix of A derived by removing all rows and columns 
with indices not in a and /3, respectively. We denote A a := A a ^ a for brevity. 
The determinant det(A a ) is called a principal minor of order |q|. 
The characteristic polynomial of J is given by 

Pj(A) = det(AI - J) = A n + Cn-iA"- 1 + • • • + Cl X + c . 

Its coefficients Cj can be computed by summing over all principal minors of 
J of order n — i: 

c i = (-l) n - i Yl det(J Q )= Yl det (" J «). ( 2 ) 

ckC {1 , . . . ,n} ctC{l,...,n} 

j ck | =tl — i | o; | — n — i 

where we set c n = 1 for completeness. Two well-known special cases are 
Co = det(— J) and c n _i = tr(— J). 

The roots of Pj are the eigenvalues of the system. A saddle-node bi- 
furcation, at which an equilibrium point splits into two, requires a single 
zero eigenvalue. The condition cq = and therefore the vanishing of the 
determinant in at least one point is a necessary condition for this bifurca- 
tion. A Hopf-Andronov bifurcation, at which an equilibrium changes into 
a limit cycle, requires a single pair of conjugate eigenvalues with zero real 
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part. Under certain conditions, the vanishing of coefficients q = for i ^ 
implies that the Hurwitz determinant of order (n— 1) vanishes, which in turn 
is a necessary condition for a conjugate pair of imaginary eigenvalues |17| . 
These observations give a direct connection of conditions for bifurcations 
and (non-)vainishing of principal minors via equations ([2]). For fully open 
systems with inflow rates K and internal and outflow reactions given by 
N, the condition that two positive equilibria match the same inflow rates is 
given by 

K = -Nv(x') = -Nv(x**) ; 

multistationarity can be excluded if — Nv(x) is injective. For mass- action 
kinetics, this requires that det(— J) > for all x > [5j, for NAC kinetics, 
the same has to hold for all its principal minors [2] 

Before introducing the interaction graph itself, we briefly recall some 
standard definitions for arbitrary directed graphs G = (V, E, 7) with vertex- 
set V, edge-set E, and edge-label function 7 : E — > M: an edge e = (u, v!) £ 
E starting in u and ending in v! is said to be incident to either vertex and 
conversely either vertex is incident to e. 

A path of length q or q-path is a sequence of edges ((«i, U2), (112, U3), . . . , [u q 
such that (ui, Uj+i) £ E for 1 < i < q. We will alternatively also denote 
it by its sequence of vertices (u%, . . . , u q ). It is a simple path if U{ ^ Uj for 
i 7^ j. It is a cycle C of length q, if u\ = u q and a simple cycle or circuit if 
it is a cycle and a simple path. A cycle (u, u) G E is a self-loop (of length 
q = l). 

A sub-graph H = (V(H), E(H)) of G, denoted H C G, is a graph with 
vertices V(JET) C V, edges f?(JET) C (V(IT) x V{H)) n J5 and edge-label 
function 7 1 K(-ff)- We define the label of the whole sub-graph H to be 

7(F) := ]J ^ e )- 

A sub-graph H is positive (negative) if 7(F) > ('y(H) < 0). Two sub- 
graphs H and ii 7 are disjoint if their vertex and edge sets are. The union, 
intersection, and difference of (sub-)graphs H, H' are defined by the cor- 
responding operations on vertex and edge-sets and restriction of 7 to the 
resulting sets, e.g., H\H' is the graph given by 

V(H\H') = V{H)\V{H') 

E(H\H') = (E(H)\E(H')) n {V{H\H') x V{H\H')) 

and edge-label function j\e(h\H')- Moreover, H + H' denotes the direct 
sum, i.e., the union of disjoint sub-graphs. 

The interaction graph is constructed by interpreting the Jacobian matrix 
as the (labelled) adjacency matrix of a directed graph. 
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Definition 3.1 (interaction graph Gj). Let J be the Jacobian matrix of a 
chemical reaction system with species Si, . . . , S n . The interaction graph 
Gi = Gi(3) = (V, E, 7) is the directed edge-labelled graph given by 



V 
E 

7(e) 



{Si, . . . , S„} 

{(S h Sj) G V x V I Jji # 0} 
Jji for e = (Si, Sj) £ E . 



Note that the topology of Gi is independent of the species-concentrations 
x, but that the edge-labels are functions in x. 

Definition 3.2 (sign-definite sub-graph). A sub-graph H C Gj is sign- 
definite if 

sign(7(H)(xo))) 

is independent of xo. 

A sub- graph is sign-definite if each of its edges is; the converse is not 
necessarily true. Under our assumptions on the network, sign-definiteness 
of an edge in Gi can be established from the stoichiometry alone. 

Lemma 3.3 (sign-definite edges). Let Gi be an interaction graph. An edge 
(S, S') £ E(Gi) is sign-definite if and only if there are no two reactions 
R, R' such that S and S' are reactants in R and reactant and product in R' , 
respectively. 

Proof. W.l.o.g., let S = S { and S' = Sj. We have 

dvi 



j((S h Sj)) = Jj,i = J2 N i- 



dxi 



We look at each individual reaction Rf. if Sj is not a reactant of Ri, then 
dvi/dxi = due to NAC kinetics. Assume Si is a reactant, then dvi/dxi > 
and Nn < 0. Then, either Sj is also a reactant, in which case Nji < and 
the contribution of Ri is negative, or Sj is a product, in which case Njj > 
and the contribution is positive. Thus, the overall sign of 7((Sj,Sj)) is 
indefinite if and only if Sj is a reactant in one and a product in another 
reaction, which can be established from the stoichiometric matrix alone. □ 



The proof of Lemma 3.3 directly suggests a method to make each edge 
sign-definite by splitting it into a positive and a negative edge. Important 
features of qualitative dynamics are preserved by appropriate choice of con- 
stants [131 E2]- However, if a sign-indefinite edge is part of a cycle, this 
procedure will result in a graph with two cycles of opposite signs. 
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Figure 1: Labelled interaction graph Gi of example network. Note that 
production reactions R4 and R$ have zero partial derivatives and are thus 
not part of the graph. 



Example 3.4. The system of Example |2 . 1 1 has Jacobian 

t — I 8A d A a dB \ 

V u dA c dB dB J 

where we denoted xa = A and Vi = Vi(xA,XB) etc. for brevity. 

The corresponding interaction graph Gj is given in Figure [T] It contains 
three simple cycles: the two self-loops C\ = (A, A) with 7(6*1) = — — 
respectively C 2 = (B,B) with 7 (C 2 ) = -cff - each of length 1, 
and the cycle C 3 = (A, B, A) with 7 (C 3 ) = ■ b^- of length 2. The sign 
of each edge is independent of xo and Gj is therefore sign-definite. 

4 Determinant Expansions and Characteristic Poly- 
nomial 

Recall that for a matrix A = (a$ j) £ M nxn , the determinant expansion is 
given by 

n 

det(A) = E(- 1 ) SiSn(7r) II a ^«' ( 3 ) 

7T 1 = 1 

where tt runs over the permutation group on {1, . . . , n} and sign(7r) is its 
sign. An equivalent expansion can be formulated in purely graph-theoretic 
terms using line-graphs. 

Definition 4.1 (line-graph [12]). Let C\, . . . , C q be a collection of disjoint 
simple cycles Cj covering each vertex of G exactly once: 

V(G) = V(d) U • • • U V{C q ) and V{d) n ^(C,) = 0, i / j . 
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Their union 

g 

L:=\JC^G 

i=l 

is called a line-graph of G. We denote by the collection of all line- 

graphs of G. 

A line-graph is also called Hamiltonian hooping |23j or nucleus [8], and 
the special definition of a subgraph in [T7] relates to the same concept. 

We again interpret A as the adjacency matrix of a labelled directed 
graph G = G(A) with vertices V = {ui, . . . ,u n } and edge label function 
j((ui, Uj)) = Aj t i. For any permutation n, the term YYi=i Ai,ir(i) * s nonzero 
if and only if all corresponding edges (u^^^Ui) exist. These edges induce a 
line-graph L of G and the product is the label of this line-graph: j(L) = 
J}Aj )7r (j). The sign of the permutation is also readily extracted from the 
graph. 

Definition 4.2 (signum of a sub-graph). Let H be any sub-graph of a 
directed graph G. Let e{H) be the number of even-length cycles in H. We 
call the number 

u(H) := {-1)< H) . 

the signum of H. 

The signum oj{L) of a line-graph L is the sign of the permutation de- 
scribed by L [12J, is independent of sign('y(H)), and is directly related to 
the signum £(-£/) proposed in [23j via oj(L) = — £(— L), where — L denotes 
the line-graph with edge-labels — 7(-). 

Lemma 4.3 (Harary [12] ^1. Let A £ M nxn with graph G = G(A). A deter- 
minant expansion of A is then given by 

det(A)= £ W (L)7(L)= £ n^HC)- (4) 

LeC(G) LeC(G) CCL 

We will focus our discussion on Gt(J), but remark that all arguments also 
hold for arbitrary square sub-matrices J Q by using the corresponding sub- 
graph G/(J a ) C G/(J) induced by V(G(J«)) = I Vi G i G a}. 
Each principal minor corresponds to a set of line-graphs of G(3 a ) C G(J), 
but each line-graph is associated to exactly one principal minor and we can 
compute each coefficient c, of the characteristic polynomial by investigating 
the line-graphs of all sub-graphs G(3 a ) with \a\ = n — i. 



Example 4.4. Consider again the system of Example 2.1 Using expansion 
via permutations, the Jacobian matrix has determinant 

,,dvi dv 2 dv 1 dv 2 dvidv 5 dv 3 dv 2 dv 3 dv 5 
det(J) = -bdg im +ac—— + a—— + c—— + — — 
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With a,b,c,d > 0, all partial derivatives are non-negative. Thus, det(J) is 
sign-definite and non-negative for all concentrations if ac > bd, whereas for 
ac < bd, the sign depends on the concentrations at which J is evaluated. 

The two line-graphs in this example are L\ = C\ U C2 and L2 = C3 with 
cycles Ci as above, and thus C{Gi) = {L\,L2}. Thus, 

det ( j )= e n^)^) 

Le{L lt L 2 } CCL 

='-»° ■(-£-£)•<-> 

+ <- 1 > , -(4!)-(^)- 

Because the edge-labels are sums of terms, the above condition for positivity 
of the determinant cannot be directly derived from the cycles alone. 




5 The Directed Species-Reaction Graph 

Analysis of qualitative dynamics via interaction graphs is considerably ham- 
pered by the fact that most networks contain sign-indefinite edges. More- 
over, edge-labels are often sums of terms containing different rate-derivatives, 
making them hard to compare independently of a species concentration x. 
These problems can all be addressed by exploiting the particular structure 
of chemical reaction systems which naturally leads to a bipartite graph with 
vertices for species and reactions. Our proposed directed species-reaction 
graph directly relates to previous definitions of bipartite graphs and to the 
interaction graph. 

Definition 5.1 (directed species-reaction graph). The directed species-reaction 
graph (DSR-graph) G = (Vs,Vr, E, A) of a chemical reaction network is a 
bipartite, directed graph with edge-label function A given by the sets of 

species vertices Vs = {Si, • • ■ , S n } 
reaction vertices Vr = {Ri, ■ ■ ■ , R r } 

edges E = Esr U Ers consisting of 
rate edges Esr = {(Si,Rj) G Vs x Vr\ dvj/dxi ^ 0} and 
stoichiometric edges Ers = {(Rj, Si) 6 Vr x Vs | Nij / 0} 

and the 

'g, X(a,b) = (S i ,R j )€E SR 

N-j, if (a, b) = (Rj,Si) G Ers 

0, else. 



edge-label function A((a, b)) = < 
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The restriction of A to E$r and Ers is denoted by X$r (a rate label) and 
A as (a stoichiometric label), respectively, such that for a sub-graph H of G, 



Xrs(H) := X\ Ers (H) = H N tj € M 

(Si^eEnsiH) 

and 

Xsr(H) := X\ Esr (H) = [] A(e) . 

eeE SR (H) 

Importantly, Xsr(H)(x) > for all x > 0. Either label remains undefined 
if the respective edge-set is empty. 



The generic DSR-graph of the network from Example 2.1 is given in 
Figure [2] This graph has an intuitive interpretation: a reaction vertex 
represents the rate of the reaction, which is positively influenced only by its 
reactant species. A change in the rate on the other hand implies a positive 
change in the rate of the products, and a negative change in the rate of the 
reactants; this is reflected by the corresponding edges. For NAC rate laws, 
the rate-edge labels are positive functions, and the label of a stoichiometric 
edge is a positive constant if the species is a product, and a negative constant 
if it is a reactant of the respective reaction. In contrast to the interaction 
graph, sub-graphs of G are therefore always sign-definite. 

Proposition 5.2 (sign- definite sub-graphs). Let H C G be any sub-graph 
of a DSR-graph G such that Ers(H) ^ 0. Then, 

sign(A(fO) = sign(X RS (H)) 

is independent o/x > 0. 

We briefly discuss some main differences of the DSR-graph to previous 
definitions of bipartite species-reaction graphs: The SR-graph of [5j uses 
undirected edges and labels them by the complex in which the species oc- 
curs. Species on the same side of a reaction form a c-pair. All possible 
orientations of the graph are considered in proofs. In the DSR-graph, these 
information are encoded explicitly in the existence and direction of edges and 
two products do not form a c-pair. Instead of c-pairs, a similar undirected 
graph in [2] labels edges by +1 or —1 to the same effect. The graph proposed 
in [T7] uses directed edges, but does not contain edges from a reaction to 
its substrate. Directed edges from a substrate to its reaction can instead 
be traversed in opposite direction, while directed edges from a reaction to a 
product cannot, which also necessitates to allow semi-cycles and paths in a 
line-graph. In this graph, mass-action kinetics is also exploited by merging 
the corresponding factor from dvj/dxi (which is the substrate molecularity 
of Si in Rj) with the stoichiometric label and using relative concentrations. 
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Rs : A 



Ra ■ ( 


) -> A 




+1 






dt'3 


( 




dA ' 


' dA . . 







-1 




—a ■ 




d 




R 2 : aB - 




< *T 





as 



i?5 : B 



+1 

Re : ^ £ 



Figure 2: Directed species-reaction graph for network of Example |2.1| 
Species vertices are given as circles, reaction vertices as rectangles. Produc- 
tion reactions R4, Rq are shown explicitly, but can be neglected for analyses. 
We again used A,B instead of xa,xb for readability. 
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The proposed DSR-graph directly relates to the interaction graph of the 
same network. Most results on injectivity and bifurcations rely on the fact 
that simple paths and simple cycles in Gj translate to simple paths and 
simple cycles in G, and we therefore emphasize this relation. 

Lemma 5.3 (relation of DSR- and interaction graph). Let G = (Vs, Vr, E, A) 
be the DSR-graph of a chemical reaction network. The interaction graph 
Gj = (V (G j) , E{G i) ,7) of that network is then found as: 

V(Gi) = V s 

E{Gi) = {(S, S') eV s xV s \3ReV R : (S, R), (R, S') e E} 

7 ((S,S'))= H(S,R))-X((R,S')). 
ReV R 

Proof. The only technicality is to see that if (Si, Sj) G E(Gi), then there is 
at least one reaction Ri such that both (Si, Ri) and (Ri, Sj) are edges in E. 
This is because 

{St, Sj) € E(Gj) < > g / < > I : $ • 

Indeed, ^ = NjjS^, which allows reconstruction of 7 from A. □ 

An edge in Gj may thus correspond to several 2-paths from Vs to Vs in 
G. Each summand of the edge-label in Gi corresponds exactly to one of the 
labels of a 2-path in G. 

Example 5.4. In the network of Example \2.1\ consider the upper- left entry 

dvi dv 3 

^1.1 = ~ a a a 

OX A OX B 

in the Jacobian matrix. The corresponding edge (A, A) in Gj corresponds 
to the 2-paths (A, R\,A) with label -a§^ and (A, R 3 , A) with label -§^. 

With J = N • (dv/dx.), the stoichiometric matrix is the incidence matrix 
describing stoichiometric edges and their labels in G, while (dv/dx) is the 
incidence matrix describing the rate edges and their labels. The (n + r) x 
(n + r) adjacency matrix of G is 

R - ( ° 

\(dv/dx) 0) ' 

All 2-paths from Vs to Vs are described by the upper-left n x n sub-matrix 
of B 2 , which is just J. By extension, an edge, a simple cycle or a line-graph 
in Gi typically correspond to several 2-paths, cycles or sub-graphs in G, 
respectively. This one-to-many mapping induces an equivalence relation on 
G. 
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Definition 5.5 (equivalence; species-cycle; species- line-graph) . Consider a 
reaction network with DSR-graph G = (Vs, Vr, E, A) and interaction graph 
Gi = (V(G I ),E(G I ),j). Let e = (S,S r ) G E(Gj) be an edge and denote by 

(e) := {((S, R), (R, S')) eExE\ReV R } 

the set of all corresponding 2-paths. Two 2-paths p,p' are Gi -equivalent if 
p,pf g (e). 

The equivalence relation is extended to cycles C = (ei, . . . , e q ) and line- 
graphs L = C\ + • • • + C q of Gi by 

(C) : = {(pi,. - . ,p q )\Pi G (e,)} 

(L) := {-Di + • • • + -D 9 | -Dj G (Q)} . 

An element in (C) or (L) is called a species-cycle or species-line-graph in 
G, respectively. Each species-line-graph is a set of simple disjoint species- 
cycles covering each species-vertex exactly once. We again denote the set of 
all species-line-graphs in G by 

C(G):= (J (L) . 

Example 5.6. Consider a Michaelis-Menten type mechanism, given by re- 
actions 

R 12 :E + S^ES 
R 3 :ES -»• E + P 

The DSR-graph for this mechanism is given in Figure [3j The edge e = 
(ES, E) G E(Gj) corresponds to the equivalence class (e) = {(ES, R2,E), (ES, R3, E)} 
in G. The cycle C = (ES, S, ES) of Gi thus corresponds to two species- 
cycles in G, using either two-path from (e) together with the two-path 
(E,R lt ES)'. 

(C) = {(ES,R 2 ,E), (E,R U ES), (ES, R 3 , E, ES)} . 

The label-function 7 of Gj can be reconstructed from the label-function 
A of G by summing over elements of the corresponding equivalence class. 
For cycles and line-graphs, this allows either a product-of-sums or a sum-of- 
products representation of the label. 

Lemma 5.7 (label- function) . Let e be an edge, C be a cycle, and L be a 
line-graph of an interaction graph Gi and let A be the edge-label function of 
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E 



R x : E + S 


->■ ES 


dvi 
~dS 




-1 


R 2 : ES^ 


E + S 








i?3 : ES — y 


E + P 




Figure 3: DSR-graph of Michaelis-Menten type mechanism of Example |5.6| 

the corresponding DSR-graph G. Then, 
7 (e) = £ A(p) 

P6(e> 

7(c) = ^ a(C) = n e A ^ = e n A ( e ) 

C"e<c) eeCpe(e) C"e(c)eec" 

7 (L) = x; a ( l/ ) = n e = e n ^ • 

i'e(i) CCLC'<=(C) L'&(L)C'CL> 

6 Non-vanishing Determinants and DSR-Graph 

We now turn our attention to the expansion of det(J) (and consequently 
det(J a )) in terms of the DSR-graph. We are particularly interested in con- 
ditions that guarantee that the determinant does not vanish for any positive 
states x > 0. 

Lemma 6.1 (signum of sub-graph). Let H C G be a sub-graph of G and 
define the signum of H as 



a(H) := (-1 



AH) 



with e(H) the number of cycles in H with even number of species-vertices. 
Then, 

u(C) = a(D) 
for any cycle C C Gj and D £ (C) . 
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Proof. A species-cycle D in G with k species- vertices corresponds to a cycle 
C in Gi of length k. Thus, u{C) = {-l) k = a{D). □ 

A determinant expansion purely in terms of a DSR-graph is now easily 
found. 

Lemma 6.2 (determinant expansion by DSR-graph). Consider a chemical 
reaction network with Jacobian J and DSR-graph G. Then, 

det(J) = ^ a ( L ) A ( L ) • 
LeC{G) 



Proof. Applying the definition of C{G) from Lemma 5.5 we reduce the 



expression to the one found for Gj in Lemma 4.3 



LeC(G) LeCid) L'e(L) L'e£(G/) 

where oj and 7 are again the signum and label function in Gi. □ 

Example 6.3. Consider the DSR-graph of Example 2.1, given in Figure [2] 
Its species-line-graphs are 



Li 
L-2 
L 3 
U 



{A,R X ,B,R 2 ,A) 

(A,R 1 ,A)U(B,R 2 ,B) 

(A,R 1 ,A)U(B,R 5 ,B) 

{A,R 3 ,A)U(B,R 2 ,B) 

(A,R 3 ,A)U(B,R 5 ,B) 



corresponding directly to the five expansion terms 

,.dvi dv 2 dv t dv 2 dvi dv Fj dv 3 dv 2 dv 3 dv 5 
d6t(J) = - bd ^^ +aC ^ A ^ +a d^ A d^ +C d^ A d^ + d^ A d^ ■ 

N V ' V V ' V v ' V v ' V v ' 

Ll L2 ^3 ^4 L5 

In contrast to the determinant expansion from Gj, the sum-of-products 
representation of labels of G yields a direct correspondence of species-line- 
graphs and expansion terms. Terms have the same partial derivatives if their 
species-line- graphs have identical substrate- reaction edges. This observation 
motivates to identify compatible line-graphs in G and determine their overall 
contribution to the expansion from their stoichiometric labels. 

Proposition 6.4 (compatibility). Let H,H' be two sub-graphs ofG. The 
relation 

H ~ H' <=^ H, H' compatible : E SR (H) = E SR (H') 
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defines an equivalence relation. We write 



[H] := {H C G\H' ~ F} 

/or #ie equivalence class of a sub-graph H . In particular, ~ partitions C{G) 
into equivalence classes in the quotient set C(G)/ ~ 

Proof. Reflexivity, symmetry, and transitivity of ~ are obvious. □ 



As an example, L\ and L2 in Example 6.3 are the only non-trivially com- 
patible species- line-graphs. The notion of compatibility suggests a strategy 
to determine if det(J) vanishes by summing over each individual compat- 
ibility class of C(G). If all classes are either non- negative or non-positive, 
the sign of the determinant is independent of the state x. 

Definition 6.5 (stoichiometric term of equivalence class). Let G be a DSR- 
graph and consider a compatibility class [L] £ C(G)/ ~. The term 

A([L]) := £ a(L')X RS (L') 
L'e[L] 

is called the stoichiometric term of [L]. It is a constant independent of x. 

The stoichiometric term can be computed from the stoichiometric matrix 
alone. 

Lemma 6.6 (computing stoichiometric terms). Fix a species-line- graph L € 
C(G). Let tj be the index of the unique reaction with substrate Sj in L. 
Define the n x n matrix by 



I, if{S h Rr.)eE SR {L) 
0, else , 



and let be the nxn stoichiometric matrix with columns not in r%, . . . ,r n 
removed. Then, 

A([L]) = det (Nl-Wi) = ±det(N L ) 
with the sign determined uniquely by W^. 

Proof. The matrix N^-W^ corresponds to a graph with only rate edges cor- 
responding to substrate-reaction pairs of [L] and all other reaction vertices 
removed; its determinant sums the contributions of all remaining species- 
line-graphs in that graph. Moreover, is simply a permutation matrix 
whose determinant is thus ±1. □ 
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Theorem 6.7 (determinant expansion by compatibility classes). Fix any 
NAC reaction system and let 3 be its Jacobian matrix and G its DSR-graph. 
Then, 

det(J) = Y, A (W) • X sd L ) • 

[L]e£(G)/~ 

The determinant is non-negative for all x > if 

A([L]) > for all [L] G C{G)/ ~ 

and positive if in addition A([L]) > for at least one compatibility class. 
Similar conditions hold for non-positivity (negativity) of the determinant. 

Proof. Because ~ is an equivalence relation on C(G), we can partition the 
sum of Lemma |6.2| by each equivalence class. Note that the rate label is 
identical and non- negative for all members of a class. If all stoichiometric 
terms are non-negative, then so is the sum. □ 

As a consequence, we can restrict our attention to finding conditions 
that establish non-negativity (or non-positivity) of the stoichiometric term 
for each compatibility class. For reasons of convenience that will become 
apparent shortly, we will focus our attention on G(—3), which is found 
by inverting signs of stoichiometric edges. With det(J) = (— l) n det(— J), 
the determinant expansion in G(3) is obviously non-zero if and only if the 
expansion in G(—3) is. We start by giving a sufficient condition to find a 
positive expansion term in open networks. 

Lemma 6.8 (existence of positive term in open networks). Consider a re- 
action network and assume that there is an inflow reaction — )■ Si and an 
outflow reaction B4 : Si — >• for each species Si, 1 < i < n. Consider the 
DSR-graph G(—3) and let L ou t be the species-line- graph 

n 

Tout := (Si, Ri, Si) , 
i=\ 

of 2-paths from each species to itself via its outflow reaction. Then, 

1- [Tout] = {Lout} 

2. A([L ut}) > 

3. A([L 0Ut ]) ■ XsR(L ut) > for all x > 

Proof. First note that Si is a substrate to Ri, so ((Si, Ri), (Ri, Si)) G E x E 
for each species Si. Thus, L ou t is indeed a species-line-graph. Since there 
is no other edge out of an outflow reaction, no other species-line-graph can 
use the same substrate-reaction pairs and thus the compatibility class has a 
single element. Each cycle in L out has exactly one species- vertex and thus 
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o~(L out ) = +1. Further, its stoichiometric edge is positive in G(— J) and 
thus the stoichiometric term of [L ut] is also. For positive concentrations, 
the outflow rate changes are positive, proving the last claim. □ 

Next, we show that species-line-graphs that either contain both the 
forward- and backward reaction of a reversible reaction always lead to a 
zero stoichiometric term for their compatibility class. This result also gives 
a retrospective justification for splitting reversible reactions. 

Lemma 6.9 (zero-contribution of reversible reaction). Let L G C(G) be a 
species-line- graph. Let Rf,Rb be the forward and backward reaction of a 
reversible reaction and assume Rf,Rb G Vr(L). Then, 

MM) = . 

Proof. The construction in this proof is illustrated in Fig. [4] (left). Denote 
by Sf,S r G Vs(L) the species with (Rf,Sf),(Rb,Sb) G Ers(L). Because 
Rf,Rb constitute one reversible reaction, the edges (Rf,Sb),(Rb,Sf) exist 
in the graph G. Construct a sub-graph V by replacing (Rf, Sf), (Rb, Sb) by 
(Rf, Sb), (Rb, Sf). If Rf, Rb are contained in one species-cycle C3 G L, they 
are now contained in two different cycles C\,Ci G L' (or vice- versa). Thus, 
V is a species-line-graph and V ~ L. Moreover, X((Rf, Sf)) = — X((Rb, Sf)) 
and \((Rb,Sb)) = —X((Rf,Sb)) and thus \(L) = X(L'). However, let rii := 
\E(Ci)\, then 

a(d)a(C 2 ) = (-l)" 1 " 1 (-IP" 1 = (-lp~ 2 = -a(C 3 ) 

and thus a(L')X(L') = —a(L)\(L). This construction gives a bijection be- 
tween species-line-graphs L and L' and thus A([L]) = 0. □ 

Lemma 6.10 (existence of positive compatible species- line-graphs). Let 
G(—J) be a DSR-graph and consider a species-line- graph L containing a cy- 
cle C. Construct a sub-graph Lc C G by replacing each 2-path ((S, R), (R, S')) 
in C by {{S,R),{R,S)). Then, 

1. a(L c )X(L c ) > for all x > 

2. L C ~C 

3. (L\C)U L c ~ L. 

Proof. The construction is illustrated in Fig. [5} With (S, R) G E(C), species 
S is a reactant of R. Then, also (R, S) G E(G) and thus Lc is a proper sub- 
graph of G with exactly one odd-length cycle Di from each species-vertex 
Si G V S {C) to itself. Thus, cr(A) = +1 and A (A) > for all x > and 
i = 1, . . . , n. These cycles use the same substrate-reaction pairs as C and 
thus C ~ Lc- They also cover each species- vertex in C exactly once, so 
L ~ (L\C) UL C . As a special case, L c = C if \Vs(C)\ = 1. □ 
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S b 



R f 






R_b 




X 








R_f 






Figure 4: Illustration of construction for Lemma 6.9 
the forward- 



A cycle containing 



and backward reaction can be split into two cycles, leading 
to a compatible line-graph with same absolute label, but opposite sign. 
Edge-labels denote stoichometric label. Dotted eges denote arbitrary paths 
through the graph, bold edges denote fixed substrate-reaction edges defining 
the compatibility class. 
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Figure 5: Illustration of compensating cycles. For a cycle Cycle 1 (C), a 
compatible compensating cycle (center, Lq) always exists, whose stoichiome- 
trics might dominate. A second compatible cycle Cycle 2 may however yield 
the same compensating cycle. Round (rectangular) nodes denote species 
(reaction) vertices. 



The sign of the contribution of a species-line-graph depends on its con- 
stituent species-cycles. Since each sub-graph of a DSR-graph is sign-definite, 
we can give simple conditions for a species-cycle to be positive or negative 
by determining the number of substrate-to-substrate 2-paths in the cycle. 

Lemma 6.11 (condition for sign of species-cycles). Consider a species-cycle 
C in a DSR-graph G{— J). Let s be the number of substrate-pairs, that is, 
of 2-paths (S,R,S') in C, such that S,S' are both a substrate of R. We call 
C a p-cycle (an n-cycle,) if a(C)X(C) > 0) (resp. <0). Then, 

C is p-cycle -<=^ s is even 
C is n-cycle <^=^ s is odd . 

Proof. We consider the four possible combinations of Vs(C) even/odd and 
s even/odd. For Vs(C) even and s odd, the number of negative 2-paths is 
also odd, thus <r(C) = +1 and A(C) < 0, as E(C) contains an odd number 
of negative stoichiometric edges. Thus, the overall contribution of C is 
negative. The other three cases follow the same reasoning. □ 

In [2] p-cycles (n-cycles) were called e-cycles (o-cycles). As a con- 
sequence of Lemma |6.10[ we can give a simple condition when a negative 
contribution to the determinant expansion is cancelled (see also [2] ) . 

Proposition 6.12 (dominating term). Let G(—J) be a DSR-graph and con- 
sider a species-line-graph L containing an n-cycle C. Let V = (L\C) U Lq- 
Then, 

a(L)X(L) + a(L')X(L') > X RS (L C ) > \X RS (C)\ . 
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We then say that Lq dominates C . We call C a stoichiometric or s-cycle, 
if \rs(L c ) = \\rs(C)\. 



The previous Lemma 6.10 and Proposition 6.12 hold the key to deter- 
mine, for all x > 0, if a determinant vanishes or not. Clearly, for a species- 
line- graph L to give a negative contribution to the stoichiometric term A([L]) 
of its compatibility class, it contains an odd number of n-cycles. Replacing 
one of these n-cycles C by Lc leads to a new compatible species-line-graph 
with positive contribution that compensates the negative. However, there 
might be a second species-line-graph in the same class, also with negative 
contribution that contains another negative cycle that leads to the same 
compensating species-line-graph. Thus, the compensating species-line-graph 
needs to dominate the sum of all these contributions. Next, we give con- 
ditions when such situation arises and provide a simple sufficient condition 
that excludes it. 

Theorem 6.13 (species-reaction intersection of species-cycles). Consider a 
species-line- graphs L in G{— J) and assume that in every compatible species- 
line-graph, each n-cycle C is dominated by Lc- Further assume that 

A([L]) < . 

Then, there are L\,Li G [L] and two non-disjoint n-cycles Ci C Lj, C 2 C L 2 
such that all paths in C% n Ci start in Vs and end in Vr. 

Proof. An illustration for this proof is given in Fig. [6j First note that the 
non-empty intersection of two cycles is always a collection of paths. Let P be 
one of the paths in the intersection of C\,Ci. If P ends in Vs, there are two 
substrate-reaction pairs involving the same substrate species. Hence, L\, L2 
cannot be compatible. If P begins in Vr, there are two different species- 
reaction vertices leading into it, one from C\, one from C2. These edges are 
contained in all compatible sub- graphs leading to non-simple cycles. The 
sub- graphs are hence not species- line-graphs. The only remaining case is a 
path from Vs to Vr which yields both unique substrate-reaction pairs and 
(potentially) the same number of reactions and species. □ 

As a corollary of the theorem, we get a simple condition of the cycle- 
structure of a DSR-graph that allows to test if the determinant of the net- 
work does not vanishes anywhere. An equivalent condition was first formu- 
lated in [5]. 

Corollary 6.14 (necessary condition for positive determinant). Consider a 
DSR-graph G(—J) of an NAC network. The determinant det(— J) is positive 
if 

1. there is at least one positive stoichiometric term 
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Figure 6: Illustration of non- feasible intersection of cycles in compatible line- 
graphs. Left: Both possible intersections ending in Vs require two substrate- 
reaction pairs with the same species and cannot occur in compatible sub- 
graphs. Right: An intersection from Vr to Vs uses a reaction vertex twice 
and leads to a zero overall contribution. Round (rectangular) nodes denote 
species (reaction) vertices. Dotted edges denote arbitrary paths through the 
graph, bold edges denote substrate-reaction pairs. 



2. every cycle C in G is either a p-cycle or dominated by Lq 

3. no two n- cycles have an intersection from Vs to Vr 

The existence of a positive term is guaranteed for open networks by 
Lemma 16.81 

7 Extensions 

A possible extensions to the idea of decomposing the determinant expan- 
sion using equivalence classes can be achieved if a partial order on the rate 
derivatives dvi/dxj can be established. Since Xsr(L) is a product of rate 
derivatives, this order induces a partial order on some compatibility classes 
L, L' such that \sr(L) > Xsr(L') for all x > 0. This order of classes implies 
consequently implies 

A([L]) > A([L']) A([L]) • Xsr(L) > A([L'}) ■ X SR (L') , 

allowing to compare two equivalence classes purely by their stoichiometric 
term. A negative contribution of [L'] might thus be compensated by a larger 
positive one of [L], independently of x. 

For mass-action kinetics, such a partial order is given for all species-line- 
graphs that use the same set of reaction vertices. First, note that the rate 
derivatives for mass-action rate laws are 

r)v ■ I 
-±(pc*) = -N j)i -v i (pc*)- — . 
dxj J x* 
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Consider now two species-line-graphs L,L' in G(— J) with Vr(L) = Vr(L'). 
Then, 



(Sj,Ri)eE(L>) 



and thus 

^Sr(L) U.(Sj,Rj)eE(L) ~ N 3,i 

^SR(L') U.(S j ,R i )eE(L>) ~ N hi 



independently of x > 



This allows to combine the results for all compatibility classes with the same 
reaction- vertices, a strategy employed in [T7]: Let Vj^ be the specific set 
of n reactions, and let G* C G(-J) be the DSR-graph with V R (G*) = V£. 
Then, all species-line-graphs in G* give a combined non-negative contribu- 
tion to the determinant if 

£ MIL])' I ft -AU>0. (5) 

[X]6£(G*)/~ \(Si,Rj)eE(L) J 



We remark that this overall term can also be computed similar to Lemma 6.6 
by replacing by a suitable nxn matrix extracted from N~ . A term ^ 
is called a fragment (of size n) in [T7] and a critical fragment, if it is negative. 
In that publication, the question of the relation between critical fragments 



and conditions on Corollary 6.14 was raised. Since the products Yl(~^i,j) 



are all positive, Corollary |6.14 gives a sufficient condition to exclude a critical 



fragment, as it establishes non- negativity for each summand in ([5]). 



8 Discussion 



The particular structure of dynamic chemical reaction network models allows 
to derive conditions to show or exclude various qualitative dynamic behavior. 
These conditions enable a first analysis and model selection independently 
of numerical values for rate constants and for all members of a large class 
of rate laws. 

In this paper, we brought forward a new definition of a bipartite species- 
reaction graph, termed DSR-graph. In contrast to previous definitions, all 
relevant features of cycles, such as feasible directions to traverse edges and 
substrate/product relations of species and reactions are directly encoded in 
the graph. Our DSR-graph contains previous definitions as special cases. 
We elucidated the direct connection of the DSR-graph to the systems' in- 
teraction graph and demonstrated how cycle features can be mapped by 
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simple equivalences of edges and 2-paths. For calculating terms of the de- 
terminant expansions of the Jacobian matrix, both graphs yield structurally 
similar formulas, but the DSR-graph allows a more fine-grained analysis of 
the terms. As a new result for bipartite graphs of chemical reaction sys- 
tems, we proposed a simple equivalence relation on the species-line-graphs 
of the DSR-graph that allows to collect comparable terms in the expansion 
and subsequently enabled simpler and more direct proofs of conditions for 
the non-vanishing of principal minors of the Jacobian matrix. We finally 
addressed a question raised in [T7] on the relation of their conditions to the 
ones developed by Craciun et al in [5]. 
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